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Abstract. A two-loop renormalization group analysis of the critical behaviour 
at an isotropic Lifshitz point is presented. Using dimensional regularization and 
minimal subtraction of poles, we obtain the expansions of the critical exponents v 
and r/, the crossover exponent ip, as well as the (related) wave- vector exponent f3 q , 
and the correction-to-scaling exponent ui to second order in eg = 8 — d. These are 
compared with the authors' recent e-expansion results [Phys. Rev. B 62 (2000) 
12338; Nucl. Phys. B 612 (2001) 340] for the general case of an m-axial Lifshitz 
point. It is shown that the expansions obtained here by a direct calculation for the 
isotropic (m = d) Lifshitz point all follow from the latter upon setting m = 8 — eg. 
This is so despite recent claims to the contrary by de Albuquerque and Leite [J. 
Phys. A 35 (2002) 1807] . 



PACS numbers: PACS: 05.20.-y, ll.10.Kk, 64.60.Ak, 64.60.Fr 

1. Introduction 

The concept of a Lifshitz point, introduced in 1975 Q, has attracted considerable 
attention during the past 25 years||. Recently there has been renewed interest in the 
critical behaviour at such points. In particular, field-theory approaches have been 
utilized to determine the dimensionality expansions of various critical indices needed 
to characterize the critical behaviour at m-axial Lifshitz points ||, [|, ||, |?], 

We assume that the order-parameter symmetry which is spontaneously broken on 
the low-temperature side of the critical line on which the Lifshitz point is located is 
O(n), and that the wave- vector instability which sets in at this point is isotropic in the 
m-dimensional subspace of M. d . A familiar model representing the universality class 
of d-dimensional systems with short-range interactions and an n-component order- 
parameter field 4>{x) at such an (m, d, n)-Lifshitz point is described by the Hamiltonian 

H = /<^{f (A,^) 2 + \ (V,0) 2 + f (V.,0) 2 + T f 2 + | M 4 } ■ (1) 

Here the position vector x = (xu,x±) has an m-dimensional parallel component cey 
and a (d— m)-dimensional perpendicular one x±. In the Landau approximation, the 
Lifshitz point is located at to = po — 0; this approximation holds arbitrarily close to 

§ For reviews and extensive lists of references, see [j^, j^, ^, |H| . The papers |i| contain lists of the 
most recent references on the subject. 
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it for values of d exceeding the upper critical dimension, which is d*(m) = 4 + to/2 
for to < 8 @. 

Interest in the model ([!]) started already in 1975 [[t| ||, [l(], |ll|. Unfortunately, 
early analyses [jl], ||, [l0| [llj based on Wilson's momentum shell recursion relations were 
rather incomplete and produced a long-standing controversy: The series expansions 
to second order in e = d*(m) — d one group of authors derived for the correlation 
exponents 77/2 and 77/4 for general values of m ^ (or for to = 1 |JTc|] ) disagreed with 
those obtained by another one [ fi"l) for the special values to = 2 and to = 6. This 
state was certainly very unsatisfactory because the class of models defined via the 
Hamiltonian ([!]) is not only interesting in its own right, but provides a natural, simple 
generalization of the standard class of n-component |</>| 4 -models (to which it reduces 
when to = 0). Moreover, besides representing universality classes, these models are of a 
prototypical nature in that they exhibit anisotropic scale invariance, offering excellent 
possibilities for investigating the question whether and under what conditions scale 
invariance might give rise to further invariances analogous to Schrodinger invariance 
or conformal symmetries fl2[ . Clearly, their thorough theoretical understanding is 
highly desirable. 

In two recent papers [| , we have been able to perform a field-theoretic two- loop 
renormalization group (RG) analysis of the model (Q) near its upper critical dimension 
d*(m) for general values of to with < to < 8. The results yield the e expansions to 
second order of all critical, crossover, wave- vector, and correction-to-scaling exponents 
about any point (to, d* (to)) on the critical line d*(m). The expansion coefficients 
involve four well-defined single integrals j<p{rn), jo- (to), j p (m), and J u (m), which we 
have been able to evaluate analytically for the special values to = 0, 2, 6, and 8, and 
for other values of to by numerical integration. An appealing feature of these results is 
that they include both isotropic cases, namely that of the usual critical point (m = 0) 
and that of the isotropic Lifshitz point (m — d = 8 — eg J. 

The requirement that our general m-dependent expressions reduce to the correct 
results for these simpler cases in which only the perpendicular or parallel parts of 
space remain provides crucial checks. These issues were briefly discussed already in 
our previous work ||, where also a direct two-loop RG calculation for the isotropic 
case to — d was announced but not described. Our aim here is to fill this gap by 
presenting details of such an analysis and a thorough discussion of the consistency 
of its results with our previous findings for general mjj We emphasize that this 
consistency is not a priori obvious because both to and d must be expanded about 8 
while maintaining the equality to = d. As we shall explicitly verify, this consistency 
holds equally well on the levels of the final e-expansion results, the counter-terms and 
the Feynman integrals. This demonstrates that the recent claim by de Albuquerque 
and Leite iflq] , who asserted that the critical exponents of the isotropic Lifshitz point 
cannot be obtained from those for general to, is unfounded. 



|| Such a calculation is analo gou s to the one for the standard |0| 4 theory (m = 0), which can be 
found in many textbooks like [M . We had originally felt that the publication of this calculation was 
unnecessary. The fact that reference [ H| got published and the reaction of many colleagues, who 
suggested that we write a clarifying paper, made us change our minds. 
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2. Two-loop results for the isotropic case m — d 

In the isotropic case m = d, only the parallel part of space remains, so that xu = x. 
It is convenient to use the rescaled variables 

<p(x) = CT(7 d/s ip(y) , y = <r~ 1/4 x . (2) 

Then the Hamiltonian implied by equation (til) for m — d becomes 



H (iso) _ J d d y 



2 (Ai/>) 2 + y (VVO + y V> + ff l^l 4 



1 CA.;.\2 , ^0 /v7/ , /A 2 , m „,.2 , .90 L;- [ 

-d/4 > _ -1/2 



(3) 



with 

,9o = o-Q a/ ' i u , A = (7(7 1/z po ■ (4) 

We use dimensional regularization to regularize the ultraviolet (uv) singularities of 
the iV-point cumulants G^* 1 = (ip ■ ■ ■ -0} cum and vertex functions T^ N > of the theory 
in d — 8 — 68 dimensions. These uv singularities can be absorbed via the re- 
parametrizations 

i> = z l Aen . .90 F d = K t8 Z g g , (5) 

and 

tq - t lp = k 4 Z t t , A - A LP = k 2 Z x A , (6) 

where txp and Alp, the values of tq and Ao at the Lifshitz point, vanish in our 
perturbative approach based on dimensional regularization. For the normalization 
constant in equation (|5|) we use the choice 

2 1 -^/ 2 r( 5 -f)r(f-2) 2 

Fd ~ f(d = 4) ■ (7) 

In order to determine the renormalization factors Z^p, Z g , Z tl and Z\, we compute 
the Laurent expansions of the vertex functions T^ 2 ', T^ 4 ', T^, and r^^ 2 at the 
Lifshitz point t = A = 0, where the latter two involve insertions of (one half of) the 
indicated operators xp 2 and (Vi/>) 2 , respectively. The free propagator with which we 
work is given by 

G(y) = [ G(q) e^y = ^ T[(d - 4)/2] y A ~ d , (5(g) = \ , (8) 
in position and momentum space, respectively. The values of the one-loop integral 

h{q) = J d d yG{y) 2 e^y (9) 
of T^ 4 '(q 1 , . . . , q 4 ) (where q = \q 1 + q 2 \, f° r example) and the two-loop integral 

h{q)= I 'd d yG(y) 3 e^y (10) 



of r( 2 ^(<7) can be read off from the Fourier transform of the generalized function \y\^, 
which is 

r(-f) 
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according to equation (2) of reference |T(|. The results are listed in Appendix A along 
with the other required two-loop integrals. We choose the renormalization factors such 
that the uv poles are minimally subtracted, obtaining 

n + 2 g 2 



ZipZt — 1 



Z$ = 1 
n + 8 3 g 

n + 2 g n + 2 



160 e 8 



n + 8 3 



0(g* 



6 e 8 



6 



9 2e 8 
n + 5 1 



5ri + 22 1 



27 24 e 8 



.9 2 + 0( 5 3 ) 



24 e 8 



9 2 + 0(g% 



and 



77 1 i n + 2 3 ff 2 , rv 3^ 



3 32 e 8 

From these results the beta and exponent functions, defined by 
(3 g (g) = nd K \ g 



and 



(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



rj p ee nd K \ \nZ p , p = ip,t,X, 

where 9 K |o denotes a derivative at fixed values of the bare variables Ao, t and go, 
follow in a straightforward fashion. They read 



and 



vt(g) -- 



-e 8 # + 



77 + 8 

~6~ 



■g z + 



4177 + 202 
1080 



g 3 + o(g 4 ) 



n + 2 q 2 , 

77 + 2 7 (77 + 2) 2 

g 



6 J 720 

7 (77 + 2) 3 2 



0( ff 3 ) 



120 



0(g 3 



The nontrivial root of j3 g (g) is 
1 

.9" 



6e 8 4177 + 202 2 , 3 



(18) 
(19) 
(20) 

(21) 

(22) 



77 + 8 5 (77 + 8) 3 8 

By evaluating the exponent functions at the fixed point g* , we arrive at the 



dimensionality expansions of the exponents 

. „, 3 (77 + 2) el _ 3 . 

1 1 n + 2 e 8 (77 + 2) (15 772 + 89 77 + 4)6! 



4 + ??*(<?*) 



77 + 8 16 



960 (77 + 8) 3 



+ 0(4), 



2 + V\(g*) 1 n + 2 e$ (77 + 2) (15 n 2 - 163 77- 2012) ef 



4 + ?7t(5*) 



77 + 8 8 



480 (77 + 8) 2 



0(e 



8^ 1 



(23) 

(24) 
(25) 
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and 



1 21 (n + 2) e 



2 40 (n + 8) 2 
Likewise we find for the correction-to-scaling exponent: 



_ dP g (g) 

UJ — 



9g 



= (41n + 202)el 
30(n + 8) 2 



The results (p3|) and (|24|) were already obtained in reference Q^, but the 
remaining ones (25)-(|2^) not. 



3. Comparison with the results for the m-axial Lifshitz point 

3.1. Critical, crossover, and correction-to-scaling exponents 

In reference [pj , which hereafter is referred to as I, we considered the critical exponents 
as functions f(m,d) of to and d, and determined their expansions in d about d*(m) 
at fixed m to second order in e = d*(m) — d. For d = m = 8 — es, we have 

e = e 8 /2 . (28) 

The series expansion coefficients of the terms of zeroth and first order in e, found in 
I, are independent of to. To check the consistency of the above results with those of 
references || ||, we therefore must merely substitute the limiting values of the latters' 



0(e 2 ) terms for m — + 8—, and use the relation (28). 

The exponents T), v, and ui given in equations (|23j), (|24|), and ( J27j ) read in the 
notation of I 

V = r)u{d,d) , v = vu(d, d) , to = uu(d, d) , (29) 

respectively. The e expansions oirjn and ^4 are given in equations (L62) and (1.64), 
and the one for uj^ follows from equations (1.66)— (1.68) , where (I.xx) denotes the 
equation (xx) of I. Note that the contribution oc lim m ^8_ j^(m)/(8 — m) = — j'i(8— ) 
which is present in the e expansion (1.66) of the anisotropy exponent 9, does not 
contribute to the ratio LO14 = uin/O at 0(e 2 ). The only term involving j</,(to) at 
0(e 2 ) is directly proportional to j^(to). Since the latter integral is of order e%, all 
contributions to the series expansion of the correction-to-scaling exponent u> to 0(e|) 
that originate from j</)(to) vanish, as they should. 

Upon substituting the values (1.86) of the integrals j<p(8— ), J«(8— ) into 
the 0(e 2 ) coefficients, the resulting series expansions to order t\ of the exponents 
( p9| ) reduce to the above results (p^), ( p4[ ) and (|27|). Likewise, the eg expansions of 
ip = (1 + ri*)/(2 + ij*) and /3 g = J^/ip that are implied by equations (1.63) — (1.65) 
for the fixed-point quantities rj*, 77*, and rj* of I agree with equations ( pB|) and (p6|). 
The remaining critical exponents, such as the specific-heat exponent a, the order- 
parameter exponent /3, and the susceptibility exponent 7, are related to rj and v via 
known scaling and hyperscaling relations. Thus, there is no need to explicitly verify 
their consistency with the results of I. 

When referring to the results of the work jl| at the bottom of p. 355 of our previous paper we 
erroneously said that these authors had obtained the expansions of the correlation exponents V12 and 
r)i4 and the correlation-length exponents V12 and U14 to second order in eg. The alert reader will have 
noticed that the cited computation was concerned merely with the isotropic Lifshitz point, so that 
only »7;4 = n and = v should have been mentioned. 
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For the case of an m-axial Lifshitz point a variety of other physically meaningful 
critical exponents can be defined, such as the exponent ?// 2 (which governs the 
momentum dependence of the inverse correlation function f ( 2 )(q||=0, <jfj_) ~ q 2 ^ 7112 at 
the Lifshitz point), the parallel correlation- length exponent Pi2, the related anisotropy 
exponent 9 = vu/v^, and the correction-to-scaling exponent u>i2 = Olou. The e 
expansions to order e 2 of all these exponents have been given in I. These latter 
exponents are not needed in the case of the isotropic Lifshitz point. The situation 
is complementary to the m = case of the standard |e/>| 4 theory, where exponents 
requiring the parallel part of space (such as rjn and vn) loose their physical significance. 
Remarks analogous to those made in I apply here: The limits d — m — > and m — > of 
those exponents that are not required in the isotropic cases m = d and m = may well 
exist, and it is conceivable that their limiting values will turn out to have significance 
for certain problems of statistical physics. We will not pursue this question further 
here. Note, however, that the limiting value of 77/2 for to — > d appears to exist and to be 
finite since its e 2 term involves the integral j^,{m) in the combination j < j > {m)/(8 — m), 
according to equation (1.61). Since jcf>{&—) = 0, the limit to — > 8— of this combination 
is nonsingular and nonzero provided that ji(8— ) 7^ (as we expect). For the same 
reason, the contribution oc j$(m)/(8 — to) to V12 in equation (1.63) remains finite in 
this limit, as do the other contributions to the 0(e 2 ) term. 



3.2. Comparison of renormalization functions 



In the foregoing subsection we have shown that the correct series expansions in powers 
of eg of the critical, crossover, and correction-to-scaling exponents of the isotropic 
Lifshitz point follow from the results of I. An analogous result holds for other universal 
quantities. This is because the required renormalization factors Z^, Z g , Z t , and Z\, 
and hence the implied beta and exponent functions (3 g (g) and rj^,, tjt, rj\ may all be 
obtained from the m-dependent results of I. 

A direct way of seeing this is to compare the Feynman integrals on which the 
analysis of the present paper is based (and which are computed in Appendix A) with 
their analogues of I. To facilitate this comparison, it is advisable to get rid of the 
choices of the factors Fd and F m>t we absorbed in g and u (the renormalized coupling 
constant of I), respectively. To achieve this goal, we must merely divide two-loop 
Feynman integrals such as I3, ^4, and I§ by the square of e$ I2 [as can be seen from 
equations ( |A.l| ) and (1.38)]. That is, the normalized integral Ii(e, 0)/[e§ /2(e)] 2 should 
be compared with the quantity {Ii(e±, 0)/[e /2(e±)] 2 } m= 8-£ 8 , e = £8 /2, where e and e± 
are unit d and d — m vectors, respectively. Utilizing equations (LB. 14), (1.86) and 

for the isotropic case, one 



(LB. 13) of I, and the corresponding results (A.l) and (A 



finds that the pole parts of these quantities agree indeed: 



h(e ± ,0) 



[^(ei)] 5 



1 

27 8 



1 



1 

12 



•0(68) 



h(e,0) 



(30) 



m=8-e8,e=e 8 /2 

Analogous results hold for the other integrals needed 
counter-terms. We have 

3? 4 



^ 2 (e)] 2 ' 

to determine the two-loop 



h(q) 

[e 8 / 2 (e)] 2 80 e 8 



•O(ejj) 



F 2 

m,e 



jo- (to) q 4 
16 to(to + 2) 



0(c°) 



1 



(31) 



m— 8 — 68,6—68/2 
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and 



-3g 2 



[e 8 J 2 (e)] 2 16 e 8 



0(eg) 



—F 2 



jp( m ) q 2 

4 m 



0(6°) 



1 



[^ 2 (e±)] 2 



(32) 



m— 8 — e 8 ,e— es/2 

where F m e = e/2( e _i_) is the counterpart of the normalization constant F^. The 
analogues of 1% and J5 we inserted on the right-hand sides of these equations can be 
read off from equations (LB. 4) and (LB. 6). Thus the pole terms found in I reduce 
in the isotropic case precisely to those obtained here. Provided we use the same 
conventions as here to fix the counter-terms, we get identical results for the Z factors 
Z^, Z g) Z t , and Z\. 

Let us briefly comment on the fact that we evaluated the integrals I 2 and I4 in I at 
a momentum q with vanishing parallel component gii and finite perpendicular one 
In the limit m — > d, the latter becomes zero-dimensional. However, this is no cause of 
concern. First of all, we should get reasonable results in this limit (as our results in 
I show) because the opposite would indicate a failure of the analytic continuation of 
the momentum-space integrals J d d ~ m q in d — m. Second, taking q = (0, q ± ) may be 
viewed as the special choice enp(iq^-x± — ax 2 ), with a — > 0+, of a test function utilized 
to determine the pole term of the corresponding distributions G(x) 2 and the Fourier 



back-transform of I4, to position space. The result G(x) 2 F n 



5{x) + 0(e°) we 



obtained in I can equally well be derived by computing the action of G 2 on other test 
functions (integrating over a finite sphere centered at the origin would do). If we set 
d = m = 8 — eg m this result, we recover the correct pole term of the isotropic case. 
Needless to say that only the (m = <i)-dimensional part 8(x\\) of the delta function 
remains for d = m, where x — xu. 

By comparing the renormalized action implied by our re-parameterizations (^|) 
and (^|) with the renormalized action of I when m = d, one can easily see how our 
two-loop results (E^)— (JT3) for these renormalization factors can be expressed in terms 
of those of I. Since to order u 2 only the terms oc u 2 /e depend on m, we can set m = 8 
when evaluating them. The desired relations then become 

Z^(g,e s ) = [Z*Z^} u=g/2 . m=S e=e8/2 + 0{g 3 ) , 



{Z 2 Z g )(g,e 8 )= Z U Z 2 Z^ 4 



=ff/2;m=8,e=e s /2 



0(<? 3 ) 



(Z+ Z t )(g, e 8 ) = [Zj, Z T } u=g/2 . m=8}e= , 8/2 + 0(g 3 ) , 



and 



{Z^ Z x )(g,e 8 ) 



Z p Z^ Z X J 2 



0(5 3 ) , 



(33) 
(34) 
(35) 

(36) 



- u—g/2;m—8,e=€8/2 

where the replacement u — > g/2 is due to the different normalization constants whose 
ratio is F m=s _ es e8 / 2 /F d = l/2 + 0(e 8 ). 

For the associated RG functions, this translates into 

g (g; e 8 ) = -€» g + WsA m = 8, e = 0) + 0( 5 3 ) . 



rn,{g) = 2i la {g/2-m = 8) + 0(g 3 ) , Vt (g) = 2r ?T ( ff /2;8) + 0( 5 3 ) , 
and 



(37) 
(38) 



r, x (g) = 2 ^{g/2; 8) - r, a (g/2; 8) + 0(g 3 ) . 



(39) 
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The reader is invited to convince himself that the results (|T^)-(|T^) for the Z factors 
and the RG functions ([l8])-(|2l]) of the isotropic Lifshitz point are indeed recovered 
from the above formulae when the m-dependent results of I are inserted into them, 
utilizing the required values (1.86) of the integrals j^m), . . . , J u {m) at m = 8. 



4. Concluding remarks 

In summary, we have shown that the counter-terms, RG functions, and the series 
expansions in powers of e$ — 8 — d of the critical, crossover, and correction-to-scaling 
exponents of the isotropic Lifshitz point can be obtained from the results of I for 
the more general case of an m- axial Lifshitz point. We have verified their validity 
by performing an independent two-loop calculation directly for the isotropic case 
d — m = 0, using well-established standard field-theory techniques. 

The fact that a proper field-theoretic analysis of the critical behaviour at m-axial 
Lifshitz points, based on the e expansion, covers also the isotropic case m = d, should 
not be too surprising: After all, the free propagator G(x) for general m ^ d goes 
smoothly over into the correct isotropic one of equation (||) as m—d — > 0. The Feynman 
graphs of the primitively divergent vertex functions are distributions constructed from 
G(x) having well-defined Laurent expansions in e. The renormalization factors are 
determined by the pole parts of the final subtractions which the primitively divergent 
vertex functions require. Clearly, their behaviour in the isotropic limit should comply 
with that of G(x). 

In closing, let us straighten out another critique by de Abuquerque and Leite 
Jj"5| . They claimed that our two-loop calculation was incomplete because we set the 
parallel component Qii of the momentum Q to zero when computing the two-loop 

integral I±{Q, K) of the graph However, the final subtractions we determined 

in this part of our computation were those associated with the |</>| 4 counter-term and 

(2) 

the vertex function r^ 2 , both of which are momentum-independent. It is true that 
h(Q 7 K) has first order poles in e that depend also on Q». These are induced by the 
divergent sub-integral >C3\ > an d are cancelled automatically by the subtractions 
provided by the one-loop counter-terms, in renormalizable theories such as the ones 
considered in I and here. Thus the convenient choice Q» = is possible. Readers who 
would like to see this re-iterated and verified explicitly in the present context may 



consult Appendix B 
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Appendix A. Feynman integrals 

Equation ( |TT| ) yields for the integrals ^ and ( |l0| ) the results 

IM-T*^ ™-}™*-* -F t £L (A.l) 
T(d - 4) e 8 
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and 



(d-6) 



h(q) = Fj 



>> r(6 - d) T(d - 4) 



4r(5-|) 2 r(|-2)r(^-6) 



p2 ^> 4 



80 e 8 



+ 



211 - 240 lng 
3200 



0(es] 



respectively. The two-loop graph ><^[j^ involves the Fcynman integral 



-74(912)93) 



q 4 \q 



9i 2 l 



■^(|g + g 3 | 



(A.2) 
(A.3) 

(A.4) 



_ t yl+e s /2 



Using a familiar trick due to Feynman we can transform the product of denominators 
into a single power by the introduction of two Feynman parameters s and t. Performing 
the integral over q then gives 

Fd r(e 8 ) 
e 8 (47r) d / 2 r0 2 

X {stq 2 12 + (l-t)q 2 -[(l-t)q 3 

By adding and subtraction the value of the integrand at q% = 0, one easily sees that 
the pole part of £1^12^3) is independent of q 3 (noting that the pre-factor of the 
integral on the right-hand side is of order e^ 1 ). We therefore set q 3 = 0, which leads 
us to the integral 



st q 12 ] 



"}' 



(A.5) 



J(d) 



[ ds s 1 -^ (1 - s) [ dtt 3 - es (l-t)- 1+£8 / 2 (1- st)- e * 
Jo Jo 



whose dimensionality expansion 

J(8 - £ 8 



(A.6) 



(A.7) 



can be determined in a straightforward fashion. Upon expanding the pre-factor of the 
integral, one arrives at the result 



J 4(gi2,g 3 ) = F c 



2 I12 



2 £ 8 r 



2e 8 



1 

£8 



1 

12 



0(e 8 ) 



(A.8) 



(2) 

As usual, the required two-loop integrals of involve just the integrals I2 and 
I4. However, we also need the integral associated with the two- loop graph (^) of 



f/y^s 2 (g,Q), where Q denotes the momentum of the inserted operator [(Vxp) 2 /2}q. 
Since the limit Q — > of this Fcynman integral is not infrared-singular and exists, 
we can consider directly the case Q = 0. The upper line with such an insertion 
corresponds in position space to the function 

n -d/2 



D{y) 



1 



128 



d-6 



6-d 



Thus the required Feyman integral can be written as 
h{q) = Jd d yG(y) 2 D(y)e^y . 



(A.9) 



(A.10) 
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It can be evaluated in a straightforward fashion to obtain 

r (n y F 2 r(7-d)r[(d-6)/2]r(rf-4) 2 (d _ 7) 

5W/ d 8r[(io-d)/2] 2 r[(d-4)/2] 2 r[(3d/2)-7] q 

' V 0(e»). (A.ll) 



16 e 8 

Utilizing the Laurent expansions of the integrals given above, one can derive the 
results (|l|)-(|l|) for the rcnormalization factors in a straightforward manner. 

Appendix B. Cancellation of momentum-dependent pole terms 

In this appendix we show explicitly for the case of the m-axial Lifshitz point that 
the momentum-dependent poles of the two-loop graph of T^ 4 ) appearing first on the 



the right-hand side of equation (B.l) are cancelled by the subtractions which the one- 
loop counter-terms provide for the divergent sub-integral inside the dashed box. It is 
sufficient to consider the problem for a one-component order parameter (n = 1). The 
lines correspond to the free propagator of the Hamiltonian (|l|) for to = po = 0; its 
explicit forms in momentum and position space can be found in equations (1.5)— (1.11). 

Zimmermann's forest formula Jl^j tells us, for any Feynman graph, into which 
subtractions for its divergent subgraphs the lower-order counter-terms translate. In 
our case, the result is particularly simple: The one-loop |</>| 4 counter-term yields the 
subtraction 



K 



><f x$j[ . (B.i) 



Here the notation 1Z [7] is used to indicate the quantity that results from a Feynman 
integral 7" [7] of a graph 7 by making the required subtractions for all its divergent 
subgraphs but not the final subtraction for 7 itself. The boxed subgraph denotes its 
singular part. This is local in position space (namely, proportional to e _1 times the 
|</>| 4 vertex). Thus the subtracted graph involves (besides a momentum- independent 
second-order pole) a momentum-dependent first-order pole of the form e _1 times the 
momentum-dependent term of order e° of the graph that results when the box is 
contracted to a point. The latter is again the one- loop graph of T^ 4 ). 



The pole terms of the difference (B.l) must be momentum-independent and local 
in position space so that they can be absorbed via the |0| 4 counter-term. From 
equation (1.38) we know that the singular part of the divergent sub-integral I2 of I4 
agrees with the pole part of ^(ei) = F rn . e e _1 , where F m ^ t is the constant (1.39) while 
ej_ denotes a unit d — m vector. Hence the difference of Feynman integrals we are 
concerned with becomes 

nmQ - K) " / £Tif {h{ " ~ K) ~ H ' J1 ■ (B - 2 ' 

We must show that the difference TZ[I^[{Q, K) — r R,[I^{Q,K)\ ( ^ _ Q is regular in e. 

This difference is given by the analogue of equation ( |B.2|) that results through the 
replacement 

1 -Of (Qf + 2qf) - 4 (0,, • Q || + q* + Q\ ) Q„ • 9|| 



Q±\ 2 + kii + On I 4 P<*± + 0±l 2 + Ivy + 0||l 4 ] [\g± + QA 2 + ° 
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(B.3) 



Critical exponents for isotropic Lifshitz points. 



11 



The subtraction has improved the uv behaviour of the integral by two powers of 
<7H , making it uv convergent at the upper critical dimension in accordance with 
Weinberg's theorem Jl8| . Hence the momentum-dependent poles — and especially the 
Q||-dependent ones — cancel, as they should, both in I^ 4 ' as well as in T^} . 
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